o 
o 






£ 



Anomalous stress relaxation in random 
macromolecular networks 



Kurt Broderix, * 
Henning Lowe, Peter Miiller and Annette Zippelius 



Institut fur Theoretische Physik, Georg- August- Universitat, D-37073 Gottingen, 
q . Germany 



,2 ' Abstract 

Within the framework of a simple Rouse-type model we present exact analytical 
results for dynamical critical behaviour on the sol side of the gelation transition. 
The stress-relaxation function is shown to exhibit a stretched-exponential long- 
time decay. The divergence of the static shear viscosity is governed by the critical 
exponent k = <j) — (3, where <j) is the (first) crossover exponent of random resistor 
networks, and (3 is the critical exponent for the gel fraction. We also derive new 
results on the behaviour of normal stress coefficients. 
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1 Introduction 



The viscoelastic properties of incipient gels have received considerable interest, 
but are still controversial. Whereas the static critical behaviour as predicted 
by percolation theory has been confirmed in experiment as well as simula- 
tion [1], the dynamic critical behaviour, and in particular stress relaxation is 
much less understood. Conceptually, the experimental procedure is depicted in 
Figs. 1 and 2: a homogeneous, time- dependent shear flow is imposed, and the 
stress relaxation is measured. Throughout the sol phase, even far away from 
the critical point, one observes an anomalous time decay of the shear-stress 
relaxation function, which follows a stretched-exponential law exp{ — (t/t*) a }. 
However the exponent values [2-7] vary over a wide range, and in some ex- 
periments even non-universal exponents are reported, depending for example 
on molecular weight or concentration of crosslinks. Similarly, the divergence 
of the static shear viscosity is still a matter of debate. Experiments yield wide 
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ranges of exponent values [8,5,9] whose origin is not clear. Even the existence 
of a single dynamic universality class has been questioned [8]. 

From a theoretical point of view, the issue is not clear either. Scaling arguments 
which are based on simple physical pictures or analogies to other random 
systems, are in conflict with each other. For example, it has been argued 
[10] that the static shear viscosity i] n of a cluster of n polymers should be 
determined by the longest relaxation time t n of the cluster. The latter has 
been identified with the time scale for diffusion of the cluster over a distance 
corresponding to its own size: t n oc R 2 r J D n . Here R n denotes the radius of 
gyration that is assumed to scale with the mass of the cluster as R n oc n l l dl 
with df the Hausdorff-Besicovich dimension. In a simple Rouse model the 
diffusion constant behaves as D n oc 1/n. The probability to find a cluster of 
size n is taken from percolation theory to yield rj oc e~ k with k = 1v — (3. Here 
e denotes the distance from the critical point and v and (3 are the exponents of 
percolation theory for the correlation length and the gel fraction, respectively. 
In another line of approach one has tried to relate the viscosity to random 
resistor networks. For example, it was suggested [11] that the conductivity of 
a random mixture of conductors and superconductors should show the same 
critical behaviour as the viscosity, s = k. 

In this Paper we discuss Rouse dynamics in an externally imposed shear flow, 
generalised to include the effects of random, permanent crosslinks. Within this 
model we have recently found [12,13] the exact result k = — (3 for the critical 
divergence of the static shear viscosity. Here is the (first) crossover exponent 
for random resistor networks, which is discussed e.g. in [14,15]. This analytical 
result is in contradiction with all previous scaling arguments, but agrees with 
that of a molecular-dynamics simulations [16]. The afore- mentioned scaling 
arguments fail, because they ignore the multi-fractal structure of percolation 
clusters - as first noted by Cates [17]. To account for multi-fractality one needs 
to introduce another fractal dimension, the spectral dimension d s , which is 
independent of /3 and v. The crossover exponent is related to d s according 
to = vdf[(2/d s ) — 1] so that our result k = — (3 is incompatible with the 
proposal k — 2v — f3. It is also incompatible with the other suggestion s = k, 
as can be seen most easily in d — 2, where duality implies s = [18]. 

Whereas for Newtonian fluids the simple shear flow of Fig. 2 gives rise only to 
a shear stress a xy , is is well known that in a non- Newtonian fluid all six com- 
ponents of the stress tensor are non-zero [19]. If the fluid is incompressible and 
isotropic, there are three independent, experimentally observable stress com- 
ponents: the shear stress a xy , the first normal stress difference a xx — a yy and 
the second normal stress difference a yy — a zz . We compute the normal stress 
coefficients \&i and ^2 and find ^2 = 0, a characteristic result for Rouse-type 
models. The first normal stress coefficient \l/i is predicted to have a much 
stronger divergence than the shear viscosity as the gelation transition is ap- 
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Fig. 1. Polymers in a shear experiment Fig. 2. Homogeneous linear shear flow 
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proached. Even though one has performed many experiments [20] on both 
the shear-rate dependence of normal stresses in entangled or (temporarily) 
crosslinked polymeric liquids and on the time dependence of the normal-stress 
response to particular shapes of shear strain, we are not aware of any exper- 
iments measuring ^i as a function of the crosslink concentration. Previous 
theoretical work [21] relates the critical divergence of \l/i at the gelation tran- 
sition to scaling properties of the relaxation-time spectrum, see also [22] for a 
recent approach in a similar spirit. 

Finally, we present a result on the long-time decay of the shear-stress relax- 
ation function. We show that this decay is described by a stretched exponen- 
tial, which is determined by the soft-mode excitations of the clusters, see [23] 
for details. 



2 Model and Observables 



We consider a fluid of N identical molecular units each consisting of L mono- 
mers. Examples are chains or rings of length L or stars with [L — 1) branches. 
Within the Rouse-type model studied in this Paper it turns out that the crit- 
ical behaviour as well as the anomalous long-time decay in the sol phase is 
independent of the internal structure of the molecular units. We therefore only 
discuss the simplest case in detail, namely molecular units which are just mo- 
nomers, corresponding to L — 1. The reader who is interested in results for 
more complex molecular units is referred to [12,13]. 

The thermal degrees of freedom are the positions Rj(t), i = 1, . . . , N, of the 
monomers, which relax in the presence of M quenched, random crosslinks, 
each connecting a pair (i e ,i' e ), e = 1, . . . , M, of monomers. Crosslinks are 
modelled as harmonic springs 
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Fig. 3. A particular crosslink realisation and the associated connectivity matrix T 
(with all coupling constants A e = 1). 

with random spring constants A e and an overall coupling strength determined 
by the length a > 0. It is advantageous to express the potential energy in 
terms of the connectivity matrix 
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A specific realisation of crosslinks is represented as a graph Q = {i e ,i' e }^L 1 
or, equivalently, by its connectivity matrix I\ A simple example is shown in 
Fig. 3. 

We consider Rouse dynamics generalised to include the effects of permanent 
random crosslinks [24-26] and an externally applied velocity field v cxt (r,£) 
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Here, Greek indices label Cartesian coordinates x,y or z. Inertial terms are 
neglected in (4), and friction with a friction constant ( occurs when the veloc- 
ity of a monomer deviates from the externally applied flow field. The cross- 
links exert a force —dU/dTLi on the monomers, in addition to a random, 
fluctuating thermal- noise force with zero mean and covariance (£f (£)£f CO) = 
2£6 a p6ij8(t — t'). Note that we have chosen units in which the inverse temper- 
ature is equal to one. 

In the sequel we will only be interested in a homogeneous linear shear flow 

v a (r,t) :=5 ax K{t)y (5) 



with a time- dependent shear rate n(t), which is sketched in Fig. 2. Given this 
shear flow, the equation of motion (4) is linear and can be solved exactly for 
each realisation of the thermal noise [13]. 

In reaction to the externally applied shear flow, the crosslinked polymer system 
exhibits stress, whose tensor components are given in terms of a force-position 
correlation [24,25] 
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Here, po denotes the density of monomers, and the initial values are to be 
taken at more and more distant times to in the past in order to ensure that 
after averaging over the thermal noise, the system has reached a steady state 
at time t. This yields for the stress tensor 



<r(t)=x(0)l+ J dt' X (t-t')K(t') 
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confer [13], where the stress relaxation function is given by 
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and the dependence of T on the crosslink realisation Q has been emphasised 
in the notation. In (8), Eq denotes the projector on the space of zero eigen- 
values of T. These zero eigenvalues correspond to translations of whole clus- 
ters. The associated eigenvectors are constant within one cluster and zero 
outside. In the example of Fig. 3 the null space is spanned by the vec- 
tors ai := 4" 1 / 2 (1, 1,1, 1,0,0) and a 2 := 2- 1 / 2 (0, 0, 0, 0, 1, 1), and one has 
Eq = a! a* + a 2 a 2 . Within the simple Rouse model the zero eigenvalues do not 
contribute to the shear relaxation because there is no force acting between 
different clusters. The only contribution to stress relaxation is due to defor- 
mations of the clusters, as can be seen from (8). The long-time decay of x wm 
be dominated by the smallest eigenvalues of T, and we expect that anomalies 
in the long-time behaviour are due to a peculiar behaviour of the density of 
eigenvalues of T for small eigenvalues. 

To complete the definition of the dynamic model, we need to specify the 
distribution of random crosslinks. Two cases will be distinguished: 



(i) Each pair of monomers is chosen with equal probability c/N, generating 
mean-field random graphs as discussed by Erdos and Renyi [27], see also 
[28] . As a function of crosslink concentration c , the system undergoes a 
percolation transition at a critical concentration c cr i t = \. For c < c cr i t 
there is no macroscopic cluster, and almost all clusters are trees. The 
average number of tree clusters per particle is given in the macroscopic 
limit by 

_ n n " 2 (2c e~ 2c ) n 
2c n! 

(ii) Crosslinks are distributed such that the cluster-size distribution follows 
a scaling law 

r n = n"V((c C rit - c)n u ) (10) 

near criticality. This case includes both mean-field random graphs and 
random bond percolation. 

For both cases (i) and (ii) we require the random spring constants A e to be 
distributed independently from each other, as well as independently from the 
crosslink positions. Moreover, the probability for very soft spring constants to 
occur shall be sufficiently small in that sufficiently high inverse moments 



P n := d\\- n p(\) (11) 



of A e are assumed to exist. 

The combined average over crosslink configurations and random spring con- 
stants will be denoted by an overbar • . Using this notation, we implicitly 
assume that the macroscopic limit iV — *■ oo, M — ► oo, M/N — > c is carried 
out, too. 



3 Static shear viscosity 



According to (7), a time-independent shear rate n(t) = K induces a time- 
independent shear stress a xy = poi]K, which is determined by the static shear 
viscosity 

,:=/d. X (.) = ^»(i^). (12) 

o 



As mentioned before, the Rouse model does not include any interactions be- 
tween different clusters so that the viscosity allows for a cluster decomposition. 
This is apparent in (12) from the block-diagonal structure of I\ We decompose 
the graph Q for a given crosslink realisation into K connected clusters {Mk\k=i 
with each cluster A4 containing N k monomers. The cluster decomposition of 
the viscosity then reads 

n{Q) = E ^(A4) • (is) 

^=1 



In order to compute the viscosity of an arbitrary connected cluster A4, we 
exploit the analogy between a random network of harmonic springs with spring 
constants A e and a random network of resistors of magnitude 1/A e , see Fig. 4. 
Since both networks are governed by linear equations, it is plausible - and can 
be shown exactly [29] - that the resistance 1Z(J\fk\i, j) between any two nodes 
(i, j) in the cluster J\f k can be expressed in terms of the pseudo-inverse of the 
connectivity matrix. This gives rise to the exact relation [13] 

k i,j€Af k 



It remains to compute the average of the viscosity over all realisations of the 
crosslinks. For convenience, the average over different crosslink realisations is 
performed in two steps: We first average over all clusters of a given size n and 
subsequently average over all cluster sizes with the appropriate cluster-size 
distribution nr n , 

oo 

V=Y, nT nVn- (15) 

n=2 



Here, rj n is the average of the viscosity over all clusters of a given size n. 
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Fig. 4. Correspondence between a network of harmonic springs and of electrical 
resistors. 
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Fig. 5. Shear viscosity (17) in units of 
(a 2 /3 as a function of c for P\ = 1. 



Fig. 6. First normal stress coefficient 
(24) in units of ((a 2 /3) 2 as a function 
of c for P 1 = P 2 = 1. 



First, we consider the ensemble (i) of mean-field random graphs. For this case, 
all clusters are almost surely trees [27], and hence the resistance between any 
two nodes (i,j) of a connected cluster is just their weighted chemical distance. 
For this reason, and due to the independence of the distribution of the A e , rj n 
is known exactly [30] 
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with Pi being defined in (11). Together with the appropriate cluster-size dis- 
tribution (9) we find [13] 
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Fig. 5 displays rj in units of Co 2 /3 as a function of c for the special case P\ = 1. 
The exact result (17) is valid for all < c < c cr i t = | and exhibits a logarithmic 
divergence as the percolation transition is approached. 

Second, we consider more general percolation ensembles which are only re- 
quired to allow for a scaling description (10) close to criticality. Here we also 
assume that all spring constants are fixed, A e = 1 for e = 1, . . . , M. The theory 
of random resistor networks [14,15] has established a scaling relation 
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for the resistance which implies [13] via Eq. (14) an analogous scaling relation 
for the viscosity (i]) n ~ n^ with <j) being the (first) crossover exponent. To- 
gether with the scaling form (10) of the cluster-size distribution, we find [13] 
for the critical exponent of the viscosity 

k = (j)-j3. (19) 



For mean- field percolation one has <p = [3 = 1, and thus k = 0, in accordance 
with (17). For 3-dimensional bond percolation high-precision simulations [31] 
for yield k w 0.71, in good agreement with recent simulations [16]. 



4 Normal stress coefficients 



The isotropic part of the stress tensor is not significant for the simple shear 
flow (5). Thus, we concentrate on the first and second normal stress differences 
Cxx—Vyy and a yy —<Jzz, respectively For a time-independent shear rate k = n(t) 
it is customary to define first and second normal stress coefficients by 

tt! := — - y , fr 2 ;= a JH - Z zz . (20) 

PO K PO K 



One deduces immediately from (7) that ^2 = 0, a characteristic result for 
Rouse-type models. In contrast, the first normal stress coefficient \l/i is non- 
zero 

1 fCa 2 \2 1 A — E \ , . 

$! = -(?—)— Tr( -±) (21) 

1 2 V 3 > N v r 2 ' y ' 



and independent of the shear rate k. Introducing the (averaged) density 



D( 7 ):=lTr{(l-£ )5( 7 -r)} (22) 



of non-zero eigenvalues of T, one gets for the crosslink average of (21) 

*i = 5 (V)/d7^. (23) 
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The inverse second moment of D was calculated in Eq. (38) of [23] for mean- 
field random graphs with the help of a replica approach. Thus, we infer the 
exact result 
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which is valid for all < c < c cr i t = \. The moments P n were defined in (11). 
The result (24) implies the critical divergence 

^~(^) 2 ^(Ccrit-Cr 3 (25) 



3 > 240 



at the gelation transition, whereas for c — ► one has 



^ = (^-) 2 fc + 0(c 2 ). (26) 



Fig. 6 displays \l/i in units of (£a 2 /3) 2 as a function of c for the special case 
Pi = P 2 = 1. 



5 Relaxation at finite frequencies 



According to (8) and (22) the averaged time-dependent shear relaxation func- 
tion x(t) is related to the density of non-zero eigenvalues of the connectivity 
matrix T by a Laplace transformation 

Xit) =p J d 7 P( 7 ) exp j- 7 ^} . (27) 



The eigenvalue density P( 7 ) has been discussed in detail in [23]: analytically 
for mean-field random graphs and numerically for finite-dimensional perco- 
lation. Here we just recall its most prominent feature for mean-field random 
graphs: D(pj) shows a Lifshits tail for small eigenvalues, as first suggested by 
Bray and Rodgers 



Din) ~ exp { - ( ^^ y 2 }, 7 10, c<±, (28) 



provided the probability density p(X) of the random spring constants A e con- 
tains no Dirac delta functions and vanishes sufficiently fast at the origin (see 
[23] for details). In the context of gelation, the small-7-behaviour (28) of D 
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gives rise to a stretched-exponential decay of the stress relaxation function 
x(t) ~ exp{ — (t/t*) 1 / 3 } for long times. In other words, it is the soft-mode 
excitations of the clusters which determine the stress relaxation at low fre- 
quencies. 



6 Outlook 



The understanding of dynamical critical behaviour at the gelation transition is 
currently far from being satisfactory. On the one hand, the experimental data 
- if existent - scatter widely and therefore do not allow for a serious check 
of theoretical predictions. On the other hand, theoretical results, as presented 
in this Paper within the framework of a simple Rouse-type model, are likely 
to be affected by shortcomings due to the neglect of certain interactions by 
the model itself. We only mention the excluded-volume interaction and the 
hydrodynamic interaction, both of which are believed to be of importance 
for stress relaxation. It is therefore one goal to incorporate effects of these 
interactions in future activities. 

Second, it is desirable to extend the semi-microscopic approach advocated here 
to the dynamics of the gel phase. This will provide information about stress 
relaxation in the gel, in particular, the critical vanishing of the static shear 
modulus at the transition. Moreover, being a cluster of macroscopic size, the 
gel admits to ask new types of questions. For instance, as initiated in [33] , one 
may examine the spatial extent of phonon-type excitations in this random 
network as a function of the excitation energy. This may reveal a localization- 
delocalization transition, a well-known phenomenon from disordered systems 
of very different kinds. 
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